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CLf 1 INTRODUCTION 



ABSTRACT 

It is shown that accurate photometric observations of a relatively high-magnification 
microlensing event (A ^> 1), occurring close to the line of sight of a gravitational 
wave (GW) source, represented by a binary star, can allow the detection of subtle 
gravitational effects. After reviewing the physical nature of such effects, it is discussed 
to what extent these phenomena can actually be due to GWs. Expressions for the 
amplitude of the phenomena and the detection probability are supplied. 
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Gravitational Waves (GWs) are predicted by General Rela- 
tivity (GR) and their existence has been indirectly proven by 
binary pulsar timing (Hulse & Taylor, 1975; Taylor, Fowler 
& McCulloch, 1979). GWs, as received on the Earth from 
any astrophysical source, produce extremely small effects 
and no GW has been detected, yet. 

The effect of GWs on some astrophysical measurable 
quantities have been the subject of studies and proposals: 
but, as a result, only upper bounds on the strength of GWs 
have been posed. Scintillation of the starlight by the GWs 
focussing of the electromagnetic radiation toward the Earth 
(Labeyrie, 1993; Bracco, 1997) have been proposed (Zipoy, 
1966) to place an upper limit on the theoretically predicted 
GW background, in analogy with the Cosmic Microwave 
Background Radiation. To detect the same stochastic back- 
ground of the GWs, the deflection of electromagnetic beams 
has been studied by several authors (Linder, 1988; Bar- 
Kana, 1996; Pyne, Gwinn & Birkinshaw, 1996) giving, again, 
only upper bounds. The same has been made in studies on 
the time-delay in lensed Quasar images (Frieman, Harari & 
Surpi, 1994). It is worth noting a proposal to find GWs from 
SuperNovae (Fakir, 1993, 1994b) and to discriminate be- 
tween different gravitational theories (Faraoni, 1996; Bracco 
& Teyssandier, 1998), observing light deflection from the 
GWs themselves. Most of these techniques rely on the ba- 
sic idea that the strength of a GW, while extremely low on 
the Earth, would be noticeably larger if detected closer to 
the source and that the deflection angle of a light beam, 
interacting with such a GW is of the order of the metric 
perturbation in the region of closest approach. However, it 
has been recently pointed out that, at least in the standard 



GR framework and up to a certain degree of approxima- 
tion, such a statement is incorrect (Bracco 1998, Damour 
& Esposito-Farese 1998, Kopeikin et al. 1999, for a brief 
discussion see also Crosta, Lattanzi & Spagna, 1999). 

In our Galaxy the most noticeable (and predictable) 
sources of GWs are binary stars (Lipunov, Postnov & Proko- 
rov, 1987), in particular, W-UMa stars (Mironovskii, 1966). 

A microlensing event occurs when, by chance, a pre- 
cise alignment between a background source and a deflecting 
mass is experienced by the observer (Paczinsky, 1986). The 
gravitational perturbation generated by a binary star very 
close to the line of sight of the event, slightly deflects light 
ray trajectories (Durrer, 1994), introducing a small distor- 
tion in the microlensing alignment. This effect translates into 
a modulation of the light amplification that is synchronous 
with the binary star. 

While we note that a similar idea has been proposed to 
describe some light variation in Quasar microlensing (Schild 
& Tompson, 1997; Larson & Schild, 2000), here we intend to 
focus our attention to the particular case of the perturbation 
by a binary star of a microlensing event, inside our Galaxy. 

The technique that will be described in the following 
sections could lead, in our opinion, to an unambiguous de- 
tection of such subtle gravitational effects from a stable 
source. Even if such a detection has to occur just during 
a microlensing alignment, the GW source can be studied 
and observed well outside the microlensing time-span. It 
is remarkable that the perturbation of the microlensing is 
not the most sensitive-to-misalignment astronomical phe- 
nomenon one can conceive. In fact, microlensing, including 
light interference effects (Deguchi & Watson, 1986; Ulmer 
& Goodman, 1995; Jaroszinsky & Paczinsky, 1995) could 
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lead to an even greater detection probability of such elusive 
phenomena. 

The role of the approximations employed to get the so- 
lutions for the light ray propagation equations is crucial, in 
order to determine to what extent the deflection angle may 
be directly ascribed to genuine GWs. Correspondingly, how 
the measurements of such small deflection angles and their 
time behaviour coud lead to new insight on gravitational 
theory and/or GW detection, remains questionable. 

Hereafter, we intend to review the contributions to the 
light deflection angle deriving from some theoretical assump- 
tions, based on different degrees of approximations, together 
with a discussion of their physical relevance. It is worth not- 
ing that, in the calculations so far published, the role of the 
first radiative term, if any, has still to be clarified. Then, 
we compute the effects of the light deflection on the mi- 
crolensing alignment and hence on the perturbation of the 
observed light -curve; we carry out a numerical example with 
some reasonable figures and, finally, we roughly estimate the 
probability to detect such an event within the framework of 
the existing campaigns to detect and follow-up microlensing 
events. 



2 DEFLECTION OF LIGHT DUE TO A 

BINARY STAR CLOSE TO THE LINE OF 
SIGHT 

The section is organized as follows: after setting some pre- 
liminary definitions, mostly of geometrical nature, we first 
give an estimate of the perturbation h determined from 
a GW source on the flat space-time Minkowskian metric. 
Then, we deal with light angle deflection calculation un- 
der different approximations: a quasi static, post-Newtonian 
case, where the speed of the bodies in the binary star are 
considered to be zero; a post-Minkowskian treatment, where 
body speeds are not negligible, using a second order light de- 
flection formalism and, finally, we report the results by Fakir 
(1994a). After expanding the deflection angle expressions up 
to any power of the impact parameter d (hence not relying 
just on the leading term of the series) we arrange them, 
with the aim to compare their relative strength and to dis- 
cuss the approximation required to establish the nature of 
the radiative component of the perturbation. 

2.1 Preliminary definitions 

Let us consider (see Fig. a double star, whose compo- 
nents have masses mi and m 2 , in circular orbit around the 
common center of mass. The two stars are separated by a 
distance p — n + T2- Special care is to be given, here, to the 
relativistic definition of centre of mass, r cm , which is well 
known in the non-interacting case; denoting by Si the i-th 
body energy, the formula is the following: 



(1) 



Formally, Eq. Q can be recovered from the classical me- 
chanics one, just replacing all the masses with their rela- 
tivistic counterparts: 



V 1 - ^/c 2 ' 



(2) 



Centre of 
mass 




To the observer 



From the source 



Figure 1. An isolated binary star, characterized by masses m\ 
and m 2 , in circular orbits with radius n and r 2 around their com- 
mon centre of mass and angular speed ui, is approached by a light 
beam with an impact parameter d. Versor k characterizes light 
propagation direction. Other geometrical quantities are indicated 
in the figure. 



where moi denote the rest masses of the two component 
stars. For interacting bodies, however, in Eq. a binding 
energy 

Gmim 2 



W b = 



2/> 



(3) 



has to be taken into account. Imitating what we have just 
done in Eq. one has to replace the mass terms in the 
classical formula with the more complicated ones (see for 
instance Landau & Lifshitz, 1971, problem 2, section 106): 

W b 

m l =m i H , i = 1, 2. (4) 

<r 

Setting the origin of our reference frame in the centre of 
mass, by definition: 



min = m 2 r 2 



(•») 



and one can define the time-dependent impact parameters 
of the two bodies as follows: 



d\ = d + n cos tot 
d 2 = d — r 2 cos tot 



(6) 



where it is assumed that the eccentricities of the orbits are 
zero (i.e. circular orbit) and the body motion is character- 
ized by an angular speed lo. One can also define the velocity 
components of the two bodies, along the approaching light 
beam direction, as Vik — vl • k, so that: 



V\k = —ujricosuit 

V2k = Wf2 COS Lot 



(7) 



where k is a versor aligned with the unperturbed light prop- 
agation toward the observer. As usual, the system total mass 
is denoted by M = mi + m 2 and its reduced mass by: 
mim 2 

u — . (8) 

mi + m 2 

The orbital period P is related both to the binary angular 
speed lo and to the binary frequency / by the usual relations: 

LO = ^ = 27T/. (9) 

Recalling Eq. the third Kepler's law 
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P 2 M 

can be cast in the other useful form: 
GM . 



2 3 
U) p 



(10) 



(11) 



Let us assume that the binary star radiates power essen- 
tially in the form of GWs at the single quadrupole frequency 
fa = 2/ (or with angular phase speed log = 2lo) with an 
intensity given by 



W 



32G 2 4 6 



(12) 



On the other hand, the energy density w carried out by a 
GW is given by 



c 3 fdhy 
16ttG Vdt) 



(13) 



where h is the dimensionless amplitude of the GW pertur- 
bation of the metric (Derouelle & Piran, 1983; Hawking & 
Israel, 1979; Weinberg, 1972). In the particular case of a 
monochromatic GW with phase angular speed log we have 



h = ho sin(uiGt) 



(14) 



and the quadratic term appearing on the right hand side of 
Eq. 1131 becomes 



(dh\ 2 hoco% 

(«) =^[l + cos(2^)] 



(15) 



After averaging it for a time span much larger than 1/fa 
one gets: 



(f) 



l2 2 
h LO G 



(16) 



and substituting it for the binary rotation angular speed 
figure 



(f) 



2h 2 Lo 2 . 



so that Eq. 113H becomes: 



C 3 h 2 LO 2 



(17) 



(18) 
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Assuming an isotropic energy distribution around the 
binary star, the energy density at a distance d from the GW 
source is: 

W 



And 2 

and, using Eq. JHJ and Eq. QSJ one finally obtains: 
, 8 GLo 2 fip 2 

From the 3rd Kepler's law 

2 G 2/3 M 2/3 
P = 



,4/3 



hence 



/(o 



(^5/3^2/3^^-2/3 



c 4 d 



(19) 



(20) 



(21) 



(22) 



A numerical approximation of Eq. I22H can be given in MKS 
units: 



ho « 4.85 ■ 10" 



d 



and in astrophysical ones: 

(P/days) 2/3 (d/AU) 
At one GW wavelength A given by 

. 27TC 7TC 
LOG LO ' 

the metric perturbation ho becomes: 
8 G"' 3 Lo^ a pM 2 ' 3 



MA) = 



7T\/5 



(23) 



(24) 



(25) 



(26) 



2.2 Calculation technique 

Hereafter, for the sake of simplicity, let us assume that the 
deflection of a light-ray due to a GW source is confined to a 
very small region, close to the point of minimum distance be- 
tween the ligth-ray and the GW source. As a consequence, 
such a perturbation is characterized by the deflection an- 
gle, since we can assume that the interaction is essentially 
concentrated at the minimum impact point. In what follows 
let us denote such an angle by a with various pedices or 
special signs, in order to distinguish the various degrees of 
approximation used to derive such a value. 

We rewrite Eq. 68 of Kopeikin & Schafer (1999) for the 
bending angle of the light ray S, under the assumption that 
the impact parameter d is negligible with respect to the dis- 
tance of the deflecting masses from the observer, obtaining: 



4G \ - rrn (1 - V ik ) 

' c 2 

i=l 



d, 



(27) 



We stress that such a result is the same obtained by Pyne 
& Birkinshaw (1993) in their Eq. 45. Furthermore, let us 
rewrite the total bending angle A</> as the sum of several 
deflection angle contributes: 



A(p = apN + Qpm + appN + Qp + . 



(28) 



where qpn refers to the Post-Newtonian formalism, that is 
to the linearized Einstein field equation, in the slow motion 
approximation or, in other words, approximating the light 
deflection angle in Eq. 1271 for Vik — » and, with a certain 
abuse of language: 



CtPM — Oi. — QpN • 



(29) 



since we are mainly interested in the new effects arising from 
the motion of the stars in the binary system under consider- 
ation. Actually, the term post-Minkowskian usually refers to 
the whole amount given by a; for a thorough discussion see 
Thorne (1987). Furthermore, a post-post-Newtonian term 
qppn is produced by second order expansion term in G, 
while ctF refers to the deflection effect claimed by Fakir 
(1994a) and criticized by Kopeikin et al. (1999) and by 
Damour & Esposito-Farese (1998). 

In order to have an idea of the order of magnitude for 
the terms in Eg. 1271 . we just note that, by replacing m and 
d by the solar values, one obtains for the leading term the 
well known deflection angle at the edge of the Sun, namely 
w 1.75". 
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The approach used here is to compute the deflection 
angle by writing the impact parameters di and a! 2 of the 
two masses in the binary star, as a function of the average 
common impact parameter d (see Eq. [SJ . The result is then 
expanded in a series of 1/d terms. The following relations 
will be used throughout this paper: 



1 



lie 
1 



l=F£ + £ =F£ +£ =F£ + 



1 =F 2£ + 3e 2 =F 4£ 3 + 5e 4 =F 6e 5 + . 



(30) 



where, for truncating the expansion up to a certain degree 
in e, some hypotheses on the smallness of e (in comparison 
with the unity) are required. 

In this way, terms of the type mir™ ±m 2 r 2 will appear. 

Finally, it is convenient to consider as a basic geometric 
configuration, the one shown in Fig.Q where the plane of the 
binary star orbit includes the straight line joining the source 
of the ray-light under scrutiny and the observer. That very 
configuration is the one which allows the maximum effect. 

After each of the following two sub-section, dedicated 
to the most relevant approximation schemes, we shall briefly 
discuss to what extent our statements have to be weakened 
for a generic geometric configuration. On the contrary, in 
Sec. 2.7 hereafter, devoted to an overview of the the results, 
we will neglect such dependencies, since we are interested in 
esteeming the order of magnitude of light ray deflection. Of 
course, these considerations should properly be taken into 
account for any specific case. 



2.3 The quasi— static field 

We can now write down the overall deflection of the light 
beam due to the two masses as the linear superposition of 
the light deflection produced by each body. We write such an 
angle as opn, to distinguish from other sources of deflections 
that we shall examine in the paper. The angle will be given 



by: 



QPN 



4G /mi m 2 \ 



which, using Eq. @, translates into: 



4G 
c 2 d 



nil 



1 + cos uit 



+ 



11)2 



1 j COS Ult 



(31) 



(32) 



Adopting the first of the expansions in Eq. I3UI . the previous 
relation can be rewritten as: 



Qpn 



4G 



(mi + m 2 )- + 



/ 2 . 2\ COS Ult 

+ {m 1 r 1 + m 2 r 2 J — — 

/ 3 3 s. COS 3 Ult 

— [m 1 r 1 — m 2 r 2 J ^ h 



(33) 



d 4 



or in full form: 

4G [mi + m 2 



QPN = 



4G r? 

" C 2 [' 



d 



'-+ 



E, 2 „ 2n\ COS 2 " Ult 

[mxrT + m 2 rt n ) ,„„ , , + 

n = l 



d 2 n +l 



(34) 



E, 2n+l 
(miri 



m 2 r 2 2 "+ 1 ) 



^ 2 n+ 2 



We just stress that the first term of the first series ex- 
pansion in brackets is not time dependent. It simply gives 
the deflection due to the whole mass of the binary, as con- 
centrated in its centre of mass. We also note that in case of a 
perfectly symmetric binary, only terms depending upon odd 
powers of d are non zero, with the d _1 term not depending 
upon the time. 

As pointed out by Kaiser & Jaffe (1997), in some oc- 
casions the simple Schwarzschild metric perturbation (the 
same effect giving the w 1.75" deflection of light at the edge 
of the Sun) can be of a similar order of magnitude. It is 
clear, however, that the latter statement can be proven only 
under some particular conditions. In fact the larger is the 
separation of a binary star the weaker is the GW strength 
and the stronger is the Schwarzschild metric perturbation. It 
has also to be pointed out that while for an oscillating mass 
the Schwarzschild perturbation goes down with d~ 2 , in the 
case of a binary source the perturbation will goes down with 
a much faster d~ 3 law. 

The first time-varying term of Eq. {MJ can be expressed 
in astrophysical units as: 



Qp N w 1.75 ■ 10 



"ii ( pi A i Ma / P2 \ 
_ 7 m @ \r q J f m \r q J 



(35) 



(d/AU) 3 

With reference to Fig. |21 we note that this effect scales 
with the cosine of the angle £. In fact the deflection above 
mentioned disappears when the impact parameter, as seen 
by the observer, during motion remains perpendicular to the 
line joining the stars of the the binary. 

2.4 A post— Minkowskian treatment 

The deflection angle of an approaching light beam by a sin- 
gle mass is given by the usual AGm/c? times d _1 in the 
perturbing mass reference frame. If the mass is moving with 
an arbitrary speed v, one further term appears, where only 
the speed component Vk along the line of sight is relevant. 
One can can think of it as an additional deflection angle qpm 
given by the linear superposition of the two moving masses 
in the binary star: 

4G /miVifc m 2 V 2 fc\ . . 

apM = -{^r + ^r) ■ (36) 

With the perturbative approach so far described, the latter 
equation can be rewritten as: 



Qpm 



4Guj cos u>t 



m\r\ 



1 + ( r i/d) cos uit 
m 2 r 2 



+ 



(r 2 /d) cos uit 



(37) 



and expanded as 
4Guj 



12 2\ COS tot 

[mir 1 + m 2 r 2 J — h 



/ 3 3\ COS LUt 

[m 1 r 1 - m 2 r 2 J — — h 



(38) 
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Figure 2. Let us define the reference plane as the one containing the centre of mass of the binary star and the line connecting the 
observer to the source. A generic configuration for the plane of the orbit may be characterized by the couple of angles (£, if)). All the basic 
calculations are carried out under the simplifying assumption that the two planes coincide. 



(mirf + rn 2 r 2 ) cos 4 ut 

+ - h . 

d 4 



The resulting relations can be rewritten in the following 
compact and exact form: 



CIPM 



4Goj 

r3 



El 2n . 2n\ COS 2 " Lot 

{mm + m 2 r 2 ) d2n + 

n=l 



/ 2n+l 2n+l\ COS 2 ' I+1 Lut 



n=l 



(m 1 r 1 



— m 2 r 2 



d 2n + l 



(39) 



In contrast to what happens for the Post-Newtonian 
contribution, the above effect depends on the cosine of the 
angle ip. This means that, while for the configuration shown 
in Fig.^both the PN and PM terms amount to the full figure 
worked out here, for a generic configuration of the orbital 
plane (see Fig. [5] for two particular sets of cases, where only 
one of the two angles is different from zero) , these two terms 
are, in general, attenuated, but they cannot simultaneously 
vanish. 



QPPN 



15ttG 2 
4c 4 d 2 



+ 



(1 + ri/dcostot) 

ml 

(1 — r 2 /dcosu;t) 2 



7 + 



(41) 



and, after substituting the expansion given in Eq. I|30|l . one 
obtains: 



15ttG 2 
4c 4 



(ml + ml) — + 



„ / 2 2 \ C°S Lot 

-2 (toiTi - m 2 r 2 J — -r= h 

i O I 2 2 . 2\ COS 2 Ojt 

+3 [mm +m 2 r 2 ) — — h 



(42) 



Contrary to what one could expect, it is remarkable that the 
first non-vanishing, time-varying term is of the same order 
as it occurs in the post-Newtonian approach, at least when 
a non-symmetric binary is considered. As before, we also 
give the complete expression for oppn: 



2.5 Post— post— Newtonian relativistic deflection 

Post-Newtonian relativistic deflection by a mass is obtained 
from the first term expansion of the Schwarzschild metric in 
the impact distance d. Of course, it is possible to go further 
and write down the deflection angle up to the d~ 2 term as in 
Epstein & Shapiro (1980), see also Ebina et al. (2000). Let 
us write this additional contribution to the deflection angle, 
considering, as before, the linear superposition of the effect 
due to the two masses in the binary: 



Qppn = 



15ttG 2 
4c 4 



(?)' 



+ 



rn 2 



(40) 



Following the same approach we used in the previous sec- 
tions, the latter equation can be rewritten as: 



QppN 



15ttG 2 
4c 4 



m 2 + m 2 
d? 



+ 



Eo / 2 2n-l 2 2n-l\ 

In { m 1 r 1 — m 2 r 2 ) 



d 2n+l 



+ 



2n 

+ 2^(2n + l)(m 1 r 1 + m 2 r 2 ) rf2n+2 



.(43) 



2.6 Radiative term cancellation in GR 

The first estimate for the light deflection angle due to GWs 
is likely to be due to Fakir (1994a). Under some specific as- 
sumptions, he obtains that, at an impact parameter d equal 
to one GW wavelength, a light ray is deflected by an angle 
qf given by: 
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(44) 



This result and its possible extensions to different values of 
d were examined by many authors (Linet & Tourrenc, 1976; 
Durrer, 1994; Fakir, 1995; Kaiser & Jaffe, 1997), with results 
only partially in agreement among themselves. 

In particular, Durrer (1994) claims qf oc h at any dis- 
tance d. Under the above assumption, replacing h in Eq. 144H 
with ho expressed in Eq. (1221 . one gets: 



OF 



12tt G 5/3 lj 5/ VM 2/3 

7f ~ 5 



(45) 



and, expressing the relationship in MKS units: 



«F 



7.63 • 10" 59 w 5/ VM 2/3 



(46) 



where qf is given in radians, while in astrophysical units 

-10 ( M /M Q )(M/M ) 2/3 



6.28 • 10" 



(P/days) 5/3 



(47) 



for deflection angle expressed in arcsec units. However, Fakir 
(1994a) points out a slightly faster decrease of the light de- 
flection with the distance d, being h oc d~ x . Kayser & Jaffe 

(1997) confirm Fakir's result for a range of d of the order of 
A but are unable to confirm Durrer's claim. In more recent 
times Fakir's results have been strongly criticized. Bracco 

(1998) was the first to point out that the dependence from 
d _1 in the deflection angle is too optimistic, although he still 
considered the d~ 3 term he found in its place, as a radiative 
one, or, in other words, intimately linked to the GW nature 
of the perturbation. 

Later, Damour & Esposito-Farese (1998) and Kopeikin 
et al. (1999), while pointing out the same result shown by 
Bracco, also claimed that the d -3 contribution is of quasi- 
static nature and has nothing to do with the radiative nature 
of GW emitted by the binary. They pointed out a radiative 
term of the deflection, but only due to the value of the metric 
perturbation at the observer's and at the source's location 
(called edge-effects), that, of course, prevents to sense GW 
fields in positions significantly closer to the very GW source. 

While the limits of these findings are going to be briefly 
discussed in the next sub-section, we want to point out that 
both Damour & Esposito-Farese (1998) and Kopeikin stress 
that their results are due to a perfect cancellation of terms 
in d _1 in the GR framework, so that any discrepancy in the 
latter can translate into a renaissance of the qf term. In 
particular, scalar GWs are likely to introduce radiative de- 
flection angles, which, under particular circumstances, could 
become comparable to the mentioned one, so that one can 
imagine to use a measure of a in order to establish the ex- 
istence of a term of the qf type (Faraoni, 1996; Bracco & 
Teyssandier, 1998; Liu & Overduin, 2000; Will, 2001). 



angles computed according to the different approximation 
schemes we examined. Thus, for the reader's benefit, we have 
collected the above results in Tab. 1. 

We note that the only non-static term in d~ power is 
the one claimed by Fakir (1994a). Neglecting such a term, 
to find an explicit time dependence, it is necessary at least 
to look at d~ 2 terms. It is worth noting that Bracco (1998), 
Damour & Esposito-Farese (1998) and Kopeikin (1999) not 
only have pointed out that Fakir's result is not reliable, but 
also that the first relevant time dependent term is in power 
of d -3 , corresponding to our first time dependent PN term. 
They correctly state that such contribution is actually a 
quasi static one and have no radiative nature. We also point 
out that Kopeikin's approach not only gives an independent 
theoretical confirmation to Damour & Esposito-Farese's re- 
sult, but also extends it to a more general setting. How- 
ever, our first time dependent PM term, clearly is of non 
quasi static nature and it is in d~ 2 power. It appears that 
such a term is missing both in Damour & Esposito-Farese 
(1998) and Kopeikin (1999) works, the first because they 
explicitely assume that the source internal motions are non- 
relativistic, so that the time-dependent external gravity field 
is quadrupolar; the second because their approximation "ac- 
counts for the static monopole, spin, and time-dependent 
quadrupole moments of an isolated system". Moreover, we 
just recall that we derived such a contribution through a 
simple series expansion in d~ n terms of an equation appear- 
ing in a later work by Kopeikin & Schafer (1999), although 
the same result is obtainable in the same way starting from a 
work by Pyne & Birkinshaw (1993), recently confirmed with 
a completely different approach by Frittelli (2003) . For a dis- 
cussion on the various approximations methods for the inves- 
tigation of GWs, the interested reader can refer to Thome 
(1987) and Zakharov (1973). 

If the PM d~ 2 term has really to be interpreted as a 
radiative or even a gravitomagnetic one is a question we do 
not intend to address here. We now define a simplified form 
for the total PM (i.e. including PN but neglecting PPN, 
hence developing the result as the first-order term in h) 
time depending light deflection angle, [A(f>], in which one 
has to take into account that: 

a) : we consider an expansion both in the GW perturbation 
h (see Eq. 12UI 1 and in the ratio A/d between the character- 
istic wavelength A (see Eq. 1251 1 and the impact parameter 
d; 

b) : just to help the reader in grasping the dependency 
structure of such an expansion, we intentionally omit all 
the trigonometric functions appearing in its time dependent 
terms; 

c) : the static (not depending on time) terms are neglected. 
As a result, one obtains: 



2.7 Are we neglecting a relevant term? 

Let us now define the average mass m = (mi+m2)/2 and the 
mass asimmetry Am = \m,2 — rni\ of the binary star. Even if 
almost all real binaries, like the one reported in our example 
(Sec. 5), are not close to be symmetric, a series expansion in 
Am/m is still possible, provided that it is not stopped at the 
first terms. Notwithstanding, the first terms can be useful 
to identify the different d~ n dependencies in the deviation 
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\a\ 

PN 
PM 

PPN X [l5nG/ (l6c 2 )] _1 
F 



Static 



mr 2 4Amr 3 
uimr 2 /c 4uAmr 3 /c ujmr 4 /c 



Static 4mAmr 



(3vr 2 w 2 ) / (v^c 2 ) X mr 2 ? 



3m 2 r 



Table 1. Weights of the various terms with time dependencies accordingly to the various development lines. An overall 4G/c 2 term is 
removed. Also note the alternate dependencies upon the system mass m and the non— symmetric term Am. 



where rjp is just a coefficient related to Fakir's claim that 
turn out to be zero in the strict GR framework. 

The validity of the further expansion in r/d, carried 
out inside each term, resides on the fact that d ~ A for the 
cases of interest here (the ones where the various terms in 
(A/d) x are comparable for various values of x) an d that r 
can be obtained through Eq, l25ll . recalling that v = up; in 
fact, it turns out that the r/d series expansion given here is 
equivalent to a series expansion in v/c. 

Such an expression can easily be re-written in the ap- 
proximation of r/d — > (small binary star approximation, 
where the star dimension r is negligible if compared to the 
impact parameter d ) as: 



lim [Ad>] 
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(49) 



which may be clearly interpreted as a series expansion 
in A/d. Furthermore, we point out a fact that, in our 
opinion, holds true in all the calculations reported here 
and in the past (by Fakir 1993, Faraoni 1996, Bracco & 
Teyssandier 1998, Bracco 1998, Damour & Esposito-Farese 
1998, Kopeikin et al. 1999). 

As one can see from Eq. 11481 . when A becomes of the 
same order than d [i.e. in the limit A/d — > 1), all the terms 
deriving from the different approaches (PM, PN ...) to solve 
Einstein field equations become comparable in strength. Just 
some numerical coefficients appear in the relative ratios, like 
the 7r factor between the d~ 2 and d~ 3 terms. That means 
that, in our opinion, further theoretical developments could 
lead to terms in A/d of order greater than 2, but still having 
a relevance in all the cases where d is of the same order of 
A. We recall again that h oc d _1 having as a consequence a 
dependency upon d with a power law steeper than a cubic 
one. Our position is just that if a certain approximation 
technique makes a radiative term disappear up to a certain 
power, it is necessary to take into account the first non- 
vanishing term (if it exists). 

That is why it could become particularly important to 
find an astrophysical case for probing the situation described 
above, that is when higher order terms become relevant. 



3 PERTURBATION ON A MICROLENSING 
EVENT 

As pointed out in a similar situation by Bracco (1997), some 
bending of light of a given angle does not necessarily trans- 
lates into effects that are simply proportional to such an 
angle but they have to be re-scaled accordingly to the in- 
volved geometry. 

Let us consider a microlensing event, where the source 
S, the lens L and the Earth as Observing point, O, are laying 
approximately on a straight line. Let us denote by Dsl the 
distance between the source and the lens, by Dlo the one 
between the lens and the observer, and by Dso = Dsl+Dlo 
the distance between the observer and the source. The im- 
pact parameter, r, is the separation between the lens and the 
straight line joining the observer and the source: it is mea- 
sured on the lens plane orthogonal to the line of sight. The 
GW is crossing transversally the line of sight at a distance 
Dgw from the Observer. Let us also assume, hereafter, that 
the GW wavelength, denoted by A, is much greater than 
both r and r_B, the impact parameter and the Einstein ra- 
dius of the microlensing event, respectively. In this way, it 
becomes negligible any differential deflection between differ- 
ent rays focussed by the lens toward the observer. 



3.1 GW source between the Lens and the 
Observer 

According to the notations explained in Fig. |3J let us define 
an auxiliary distance p on the plane defined by the observer 
and perpendicular to the microlensing alignment line. Since 
the deflection angle is very small, we write 



p = Acf) D GW 

and, using a = p/Dso and Ar = a Dsl, one gets 



A(j> 



Dgw Dsl 
Dso 



(50) 



(51) 



where Ar is the variation of the impact parameter r, due to 
the GW. Of course such a parameter oscillates, in the case 
of a monocromatic GW, between r — Ar and r + Ar. We 
just note that when Dgw = 0, Ar vanishes, so that any GW 
source close to the observer does not produce any effect. The 
maximum effect holds when Dgw = Dlo, that is when the 
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Source 



Dso 



Lens 



Binary 



GW 




DSL 



Dlo 



Figure 3. A binary star as a source of GWs is between the lens 
and the observer during a micorlensing event. GWs will pcrtur- 
bate the alignment leading to some noticeable signature in the 
observed lightcurve. In order to make the figure clearer the de- 
flection due to the lens is not shown. 




Observer 



GW 



Dlo 



Figure 4. As in Fig.|I] but with the GWs source between the 
source and the lens. Again, the deflection due to the lens is not 
shown. 



GW is generated in the neighbourhoods of the lens. In such 
a case, the following equation holds true: 



Ar = A<j> 



D lo D S l 



(52) 



Dso 

When D = D LO = D S l we obtain Ar = A0-D/2, so 
the maximum optical lever is equal to one fourth of the 
whole distance between the source and the observer. It is 
worth noting that, in the case of Galactic measurements, this 
condition poses a somewhat upper limit on the maximum 
lever, given by roughly one half of the distance of the Sun 
from the Galactic center. 



3.2 GW source between the Source and the Lens 

With reference to Fig.^Jit is useful, in this case, to introduce 
the displacement q, given by 



q w A<?!> (D so - D GW ) 



(53) 



as in the previous subsection, let us define the angle (5 
q/Dso and the variation Ar sa /3Dlo, thus obtaining 



Ar = A<j> 



-Dso — Dgw 
Dso 



-Dlo- 



(54) 



It is easy to see that the above relation is close to Eq. 1511 . 
The behaviour of the optical lever is similar, provided that 
one replaces the distance of the GW from the Observer with 
the same distance, but measured from the source. It is also 
evident that Eq. I54H and Eq. 151H give the same value when 



the GW is located close to the lens. The same considerations 
described in Section 3.1 hold true here. 

With reference to Fig. |S] one can introduce a lever 
length Zgw, such that Ar = A<^£gw- The behaviour of £gw, 
measured on the microlensing straight line connecting the 
source and the observer is the following one: starting from 
the zero value at the source, it linearly grows since it reaches 
its maximum at the lens position, then it linearly decreases 
to zero, approaching the observer. Provided that a certain 
sensitivity to Ar is accomplished, a bi-conic volume of the 
Galaxy is probed in searching for GWs larger than a given 
threshold (see the upper right insert in Fig. |SJ . 



4 PHOTOMETRIC EFFECTS ON 
MICROLENSING 

Following Paczinsky (1986) we use the Einstein radius te 
defined as 



4GM L -D L o-D L 



D s 



(55) 



where Ml is the mass of the lensing object. Denoting by u — 
r/tE, a dimensionless impact parameter, the amplification 
factor A(u) is given by: 



A(u) = 



u 2 + 2 



(56) 



(57) 



uVu 2 + 4 

The GW within the line of sight of a microlensing event in- 
troduces a tiny perturbation in u and a corresponding mod- 
ulation of the amplification, which can be esteemed to be: 

OA 8 

~du ~ ~ v?(v? +4)Vu 2 +4 ' 

We point out that in the cases of interest here, u < 1 holds 
true, hence the previous relations simplify into A ~ u _1 and 
dA/du « ~A 2 . 

Due to the perturbation, Ar the dimensionless impact 
parameter is perturbed by an amount An = Ar/re and the 
magnification A will exhibit, in the approximation A ^> 1: 



AA 



-A z Au. 



(58) 



Such a perturbation leads to a hopefully measurable bright- 
ness variation A7 of the observed flux / (see also Fig. 0. 
Let us define Am as the maximum photometric magnitude 
difference between several measurements affected by an in- 
tensity I ± AI, i.e.: 

Am=|log(l + 2^) , (59) 
so that one can approximately write: 

By using Eq. 1581 and the concept of effective length Igw 
defined in the previous section, one can write 

5 A(j>A l G w 



Am : 



In 10 



(61) 



Taking into account Eq. H22H . Eq. H44H and Eq. 16H it is 

possible to estimate the maximum projected distance d max , 
around which a binary star produces a perturbation, leading 
to a given photometric amplitude Am: 
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Source Lens Observer 

Figure 5. The ratio between the effect on the impact parameter, 
Ar, and the deflection angle generated by the GW perturbation, 
A</>, is here defined as a lever length £qw which is maximum near 
the lens and drops linearly to zero both toward the source and 
the observer. For a hypothetical event, where the lensing object 
is located in the bulge of our Galaxy (see upper-right inset), the 
searching for a GW signature as described in the text is equivalent 
to probing for GW sources within a biconical volume. 




o I , i i i i i i i i i i i i i i i i i i i i , , 

-6-4-2 2 4 6 

Time [arbitrary units] 



Figure 6. An example of a (strong) signature of GWs generated 
by a binary star located close to the line of sight of the micro- 
lensing event. The light from the binary source (that would be 
likely to appear as an eclipsing binary) is not included in this 
plot. 



nVEn 2 g 5/3 Aigw ^Vff 

max ~ In 10 ' c* ' r E ' Am ' (b > 

where we have grouped on the right side of Eq. iUi2l the nu- 
merical coefficients, the physical constants, the microlens- 
ing parameters and the binary star's one into four different 
fractions. The linear dependence upon Igw we have found 
deserves, however, some specific comments. We just note, in 
fact, that Eq. (1551 can be rewritten as: 

te = \J r s l^ , (63) 

where rs is the Schwarzschild radius of the lens and Z M is the 
effective length of the microlensing setup, denned similarly 
to law ■ While Igw and Z M are a priori fully uncorrelated, it is 
clear that by selection effect, the largest Am can be obtained 
when Igw ~ Ip. Under this condition, the dependence in 
Eq. 1621 from such a characteristic length becomes weaker, 
i.e. a square root. Also the dependence upon the mass of 
the lensing object becomes an inverse square root. 



5 AN EXAMPLE WITH A W-UMA BINARY 

Following Mironovski (1966) we take a typical W-UMa bi- 
nary as the average of the ones listed by Kopal (1959), ob- 
taining mi = 1.46M Q , m 2 = 0.78M Q and P = 0.3 d , leading 
to w ~ 2.42 • 10 _4 s _1 . This kind of double star has a typ- 
ical size of ~ 10 _2 AU and a radial velocity v such that 
v/c « 10~ 4 . A GW wavelength A « 26AU is obtained. We 
just note that such a figure is at least an order of magnitude 
larger than a typical t-e for microlensing in the Galaxy, so 
that the approximations used in the calculation reported in 
this work appear reasonable. 

At a distance d — A, an ho ~ 1.33 ■ 10~ 15 is obtained. 
Assuming that this W-UMa is in the bulge of our Galaxy, 
the same perturbation on the Earth will be lowered to a 
mere h » 1.68 ■ 10~ 23 . 

According to Eq. 1491 1 . one estimates a light deflection of 
A<j> « 9.5 • 10" 16 , equivalent to A^ « 0.20 nano-arcseconds 



(mainly from the PM term). Here and in the following con- 
siderations, we assume that what we called Fakir's contribu- 
tion vanishes: rjF = 0; nonetheless, we note that our results 
could increase by a factor ~ 20 in case that ijf = 1. 

Let us suppose that a microlensing event occurs and the 
corresponding magnification is relatively high (A = 100); let 
us also suppose that the lens determines a value r E = 0.1AU 
and the average W-UMa under investigation is located near 
the Galaxy center, at roughly lOkpc from the Sun. With 
the further assumption that the source is situated lOkpc 
far away from the lensing mass, a Igw = 5kpc is ob- 
tained. Combining these figures in Eq. 1611 one can estimate 
Am « 2.2 • 10 -3 . Nevertheless, one should note that a pre- 
cision of the order of 1/100 of magnitude should be enough 
to detect such a GW event, provided that one takes the av- 
erage of several GW-induced oscillations around the time 
corresponding to the magnification peak. 

The assumption we made about te is roughly one order 
of magnitude smaller than a typical one. While we note that 
this is just a factor three times smaller than the typical mass 
for microlensing, we point out that similar results can be 
obtained with a smaller impact factor d. In particular, the 
PN term, growing with d~ 3 , can make a similar perturbation 
on the microlensing event when, with a more typical r_g = 
1AU, the impact factor drops toads; 5.8AU, while the PM 
term, much larger by a factor 7r at the one GW wavelength 
distance, behaving just as a d~ 2 power, becomes of the same 
magnitude at d « 5.1AU. 

All the calculations reported in this paper are carried 
out under the assumption that the rays focussed by the mi- 
crolensing are approximately subject to the same gravita- 
tional deflection. Even if this condition fits very well the 
astrophysical cases described here, of course it is no longer 
true when the binary itself is responsible for microlensing. 
It is easy to realize that in this latter case, the effect should 
be much smaller, because the focussing deflection would 
slightly change and the relative variation in the amplifica- 
tion A should be of the same order of magnitude as the ratio 
between the deflection due to the time-dependent terms and 
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foreground star 



binary star 




Figure 7. In a reference frame with the origin fixed with the 
mass causing microlensing, appearing in the centre of the figure, 
the probability that simultaneously a suitable binary is within a 
circle of diameter <i max and a foreground star is within another 
of diameter cfi ens ing is computed in the text. Moreover, the case 
corresponding to independent objects is first evaluated, and then 
the one in which the mass responsible for the microlensing is a 
third companion of the binary is considered. 



the deflection caused by the static contributions. Neverthe- 
less, since such a possibility was beyond the scope of our 
paper, we have not investigated it in more detail, so that 
we cannot exclude that further interesting results could be 
obtained along this way. 



6 DETECTION PROBABILITY 

From an inspection of the MACHO project microlensing 
alert Web page 1 we found an average detection rate of 
events with A > 8 of the order of w 5 events per year, 
with an average light amplification A ~ 20; we dropped 
from our estimate caustic crossing events, whose treatment 
is beyond the scope of this paper. In the MACHO program 
N, w 4.3 ■ 10 s stars are photometrically observed (Alcock 
et al., 1995) hence we can define an Event Detection Rate 
(hereafter denoted by EDR) r/, expressed in detected events 
per star, per year, for such a type of high-magnification 
events, as 



I) : 



N(A > 8) r 
Because A 



1.2 • 10" 



(64) 



~ one can easily use the estimation of 
r\ to evaluate the EDR corresponding to the case in which 
another star (the lens) falls within the line of sight of a given 
star in the bulge (the source) at a certain distance diensing, 
given by 



file 



^ « 0.05AU . 
A 



(65) 



In fact, such an EDR for a diameter of the order of 
dmax will be obtained by simply scaling r\ by the ratio of the 
cross-sections defined by the two distances under study (see 
Fig- 0, leading to: 

1 http : //darkstar . astro . Washington . edu 



lcnsing 



(66) 



The EDR to have both a microlensing event (like the ones 
considered here) and, simultaneously, a third object within 
a much larger distance d max can be written as the product 
of the two terms: 



1Z « 7?7 



(In 



lcnsing 



(67) 



where 7 is a coefficient that takes into account that the 
two events should simultaneously occur. Indeed, a given 
star spends a fraction 7*7 of the time undergoing high- 
magnification events, where 7 is the average duration of a 
high-magnification event, expressed in years. A typical value 
for 7 is of the order of 10 -2 . 

Let us estimate ri max for a typical Galactic event, where 
the lens is located in the bulge: in this case we assume 
law ~ 5kpc and te ~ 1AU. In order to achieve the high- 
est possible photometric accuracy, it is to be recalled that 
one can average a few GW periods in the time span cov- 
ered during the high-magnification event. Assuming a final 



error of the order of Ar 



10~ J mags (Frandsen, 1993; 



Gilliland & Brown, 1992) one can obtain, using Eq. I|61^ 
and Eq. (67), together with the average W-UMa as from 
the previous section, and the above mentioned estimates, 
a figure for d ma x ~ 6AU. We just note this is significantly 
smaller than A so that, in such a regime, some still unveiled 
terms in higher power on A/d can even dominate. 

Finally, one have to further select only the cases where 
the second star is of W-UMa type. The population density of 
W-UMa is almost constant all the way to the Galactic bulge, 
amounting to roughly pw-UMa ~ 1/280 (Rucinski, 1994, 
1997). Then, one can combine all these EDR parameters 
into a single relationship, giving the average GW-detection 
time interval r ~ ~RT X between two successive events: 



H 2 

"'lcnsing 



7?/ 



-(/- 



<. /OW-UMa. 



4.3 x 10 4 yrs. 



(68) 



This figure could make any reasonable search for such 
event, truly hopeless. It has been obtained assuming that the 
W-UMa and the lensing star position are fully uncorrelated 
(which translates into the square power dependency upon 
r\ and to the appearance of 7). However, one should also 
consider the possibility of having a triple star system com- 
posed by a W-UMa and a third companion at a much larger 
distance, responsible for the microlensing. Because d max is 
so small, we found that the triple star case, in spite of its 
unlikeness, has an even higher EDR than the simpler case 
previously discussed. 

In the triple stars case, in fact, the EDR of a microlens- 
ing event on the companion of a W-UMa, with high magni- 
fication rate, is simply given by: 



n' 



■ 7?pw-UM a 773 



(69) 



where 773 is the fraction of W-UMa exhibiting a third com- 
panion at a projected distance of the order of magnitude 
similar to ci max . Such a number turns out to be of the or- 
der of 773 « 0.2 (Herczeg, 1988; Tokovinin, 1997), leading to 
a typical time interval between two different GW potential 
detections (or at least of the gravitationally induced effects 
discussed in the text) of the order of: 
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r = w « 273 yrs , (70) 

that is a factor roughly two hundred times smaller than in 
the case of free floating W-UMa and microlensing star. We 
note that in the triple star case the condition Lgw = Ifi is 
automatically fulfilled. 

Care should be taken into account in considering such 
figures. Most of the parameters involved are rough estimates 
and several uncertainties of the order of at least a factor two 
can be considered. This leads, in our opinion, to a significant 
diminishing in the order of magnitude for t . 

One should also consider that such a time will be low- 
ered by future improved photometric capabilities, and that 
it has been devised for an average W-UMa: in principle one 
cannot exclude that, by chance, some stronger event may 
occur. Moreover, we point out that the simultaneity condi- 
tion discussed in the case of incorrelation between the W- 
UMa and the lensing star, is well verified for the more opti- 
mistic assumptions. Notwithstanding, it can introduce some 
further augmentation of r for the more pessimistic calcula- 
tions. The relatively low time-scale, however, also suggests 
that there is a small chance that a gravitational effect signa- 
ture (to be ascribed to a genuine GW effect or mainly to a 
Post-Minkowskian contribution is here beyond our purpose) 
could even be hidden in some of the microlensing photomet- 
ric runs already collected from the ground. 



7 CONCLUSIONS 

In the writers' opinion, a critical review of the light deflec- 
tion contributions due to the gravitational field of a binary 
star shows that terms directly linked to a GW (or radia- 
tive ones) cannot be excluded by current approximations, 
available in the literature. When the impact parameter is of 
the same order of magnitude of the GW's wavelength, such 
terms could be as important as the leading, non radiative 
ones. Moreover, a non quasi-static term depending upon the 
inverse square of the impact parameter is found. Such de- 
flection angles are so small that one would be tempted to 
leave them in the academic rather than the experimental 
realm. However, we have shown that such deflection, when 
occourring close to the line of sight of a micro-lensing event, 
can give rise to detectable effects. The probability of such an 
alignment is, however, extremely small, except for the case 
of microlensing by a third wide companion of a close binary 
star. In that case, from the probability to detect such an 
event, one can reasonably expect that those effects could be 
observable in the near future, or even that one or a few of 
them are hidden in the existing literature data. 

We avoided to speculate, in the previous sections, on 
some extreme cases where the gravitational effects could 
be very large. One could, in principle, conceive a geometry 
where a small lens (for instance with a mass of the order of 
Jupiter) produces a highly amplified (A w 100) microlensing 
event for a duration of several days, with a massive binary 
like /x-Sco, close to the line of sight. In this case, the bi- 
nary star could form with the microlensing event an angle 
sufficient to photometrically study the binary star well sep- 
arated from the microlensing event, for a reasonable number 
of binary periods. Even if such an event appears to be very 
unlikely, it could determine many insights on gravitational 



and GW studies, including the possibility to probe gravita- 
tion theories that do not lead to a perfect cancellation of the 
term claimed by Fakir. 

On the other hand, we believe remarkable that ordi- 
nary events, where a gravitational or GW signature can be 
clearly identified, have a significant chance to pop-up in the 
current surveys. Post-microlensing studies should permit to 
eventually confirm the GW nature of the waviness detected 
during microlensing. Such approaches include spectroscopic 
and photometric studies (to ensure in detail the nature of 
the binary star responsible of the GWs) and, in future, as- 
trometric studies carried out with some high angular reso- 
lution tool (speckle or adaptive optics) in order to establish 
the exact geometry of the three objects involved (source, 
lens and binary stars) at the moment of the microlensing 
peak magnification. 

A positive detection, such as the one described here, can 
be seen as a sort of astrometric detection of a GW effect us- 
ing a Galaxy-sized telescope whose objective is made by the 
gravitational lens responsible of the microlensing (Labeyrie, 
1994). 
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